Mathematical Models of Physical 

Systems 

Why mathematical models of physical systems 
needed? 

Design of engineering systems by trying and error 
versus design by using mathematical models. 

Physical laws such as Newton’s second law of 
motion is a mathematical model. 

Mathematical model gives the mathematical 
relationships relating the output of a system to its 
input. 



Mathematical Models of Physical Systems 

(continued) 

• Control systems give desired output by controlling the input. 
Therefore control systems and mathematical modeling are 
inter-linked. 



• Mathematical models of control systems are ordinary 
differential equations 

• To deal with linear o.d.e, Laplace Transform is an efficient 
approach. 



What is Laplace Transform? 


• The Laplace transform of a function f(t) is defined as: 

oo 

F {s) = f\f (t)\ = \ f (t)e “dt 

0 

• The inverse Laplace transform is defined as: 
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| F(s)e s/ ds 
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where j = V-T and the value of ct is 
the singularities of F(s). 



determined by 



Why Is Laplace Transform Useful ? 


Model a linear time-invariant analog system as an 
algebraic system (transfer function). 

In control theory, Laplace transform converts linear 
differential equations into algebraic equations. 

This is much easier to manipulate and analyze. 



An Example 

• The Laplace transform of e can be obtained by: 



Linearity property 


These are useful properties: 

l[kf{t)\ = kt[f(t)\ = kF(s) 

e[fi(t) +fm= arm + arm = fm + F its ) 


Name 

Time function f(t) 

Laplace Transform 

Unit Impulse 

8(t) 

1 

Unit Step 

u(t) 

1/s 

Unit ramp 

t 

1/s 2 

nth-Order ramp 

t n 

n!/s n+1 

Exponential 

e-at 

l/(s+a) 

nth-Order exponential 

t 11 e~ at 

n !/( s+a} 1+1 

Sine 

sin {bt) 

b/(s 2 +b- ? ) 

Cosine 

cos {bt) 

s/(s 2 +b /) 

Damped sine 

e~ at sm(bt) 

b/((s+af+b?) 

Damped cosine 

e at cos{bt) 

(s+a)/((s+af +b 2 ) 

Diverging sine 

ts\n(bt) 

2bs/(s 2 +b 2 ) 2 

Diverging cosine 

tcos(bt) 

(s 2 -b 2 ) /(s 2 +V f 









































TABLE 2.2 Laplace transform theorems 


Item no. Theorem 

Name 


table 02 02 



Find the Laplace transform of f(t)=5u(t)+3e 


• Solution: 



Another Example : 


d 2 


dt 


y . + 12^ 

dt 


2 



Partial Fraction Expansion 

(Case 1 : Roots of Denominator are Real and distinct) 

Find the inverse Laplace transform of 

F(s )= 5/ (t s 2 + 3s +2 ). 

Solution: 



Exercise: Do example 2.3 of the textbook 

Laplace Transform solution of a differential equation 






Note that roots are -p t 
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Example 
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s + Pl ) r F(s ) 
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Case 2: continue of the example F(s )=(2s +3)/(s 3 +2s 2 +s ). 



Case 3: Roots of the Denominator are Complex. 

Example: F(s )=1 0/( f +4s 2 +9s+ 10) 


F(s) 
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fit) 

fit) 


2e 


It 


2 e 


2 1 


+ 2.23 c' cos( 2 t - 153 .4° ) 

- 2 e~ r cos( 2t) + c r sin( 2/) 


For the Tutorial next week, use Matlab to 
find the Inverse Laplace Transform for the 
case where the roots of the denominator are 
complex roots. 
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The nth - order differential equation : 

d n j 77—1 i i m 

c a c dc 7 d r 7 

h Cl i — + . . . + Cl, h Ctr,C — u h u 

dt" dt dt dt m 

Laplace transform : 
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The characteristic equation is defined as : 

s" +a ,s n ~ l +... + a,s + a fl =0 

77—1 10 

The roots of this equation are called poles of the system 


The roots of the equation : 

bs m +b ,s m ~' +... + b,s + b t =0 

777 777-1 1 0 

are called zeros of the system. 




Transfer Function for a differential equation 


Exercise: Do examples 2.4 and 2.5 on the Textbook. 



Models of Electrical Circuits 



• Inductance circuit: 






Ls if / ( 0 ) 0 




Models of Electrical Circuits 


• Capacitance circuit: 



£ (output) 
l(input) 


= Transfer function 


(r)d r + v(0) 

O | v(0) 

5 s 

if v(0) = 0 



Transfer Functions for Electrical Networks 


TABLE 2 J Voltage-current, voltage-charge* and impedance nela p$ Tor capacitors, resistors^ and inductors 


Impt-damc Admit tanct 

Component Voltage-current Current -voltage Voltage-charge Z(s)= V(s)fl(s) y(s)=/(s)/V(s} 




Kirchhoff’ s voltage law: 

The algebraic sum of voltages around any closed loop 
an electrical circuit is zero. 


Kirchhoff’ s current law: 

The algebraic sum of currents into any junction in an 
electrical circuit is zero. 



Electrical network (Example 2.6) 




Original network 


Ls R 




V c (s) 


Transfer Function network 


Laplace Transform network 






Exercise: Do example 2.10 of the textbook 


Excluding Operational Amplifiers 
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Models of Mechanical Systems 

Mechanical translational systems. 

Newton’s second law: 








Device with friction (shock absorber): 


m = bi 



dt 


dxi{l ) 

dt 


B is damping coefficient. 



Translational system to be defined is a spring (Hooke 5 
law): 



K is spring coefficient 





dree -velocity, force-displacement, and impedance translational relationships 
for springs, viscous dampers, and mass 


Component 


Force- velocity 


F ore e-tl is phi cement 


Impedence 

Zm(s) = F(s)/X(s) 


Sprin 




K ffi V ( r ) dT 


Kx(t ) 


Viscous damper 


Mass 


<Px(t) 


table 02 04 



Model of a mass-spring-damper system: 




Note that linear physical systems are modeled by linear 
differential equations for which linear components can be 
added together. See example of a mass-spring-damper 
system. 



Transfer Functions for Mechanical Systems 




(«) 


KX(s) 
f v sX(s) 
Ms 2 X(s ) 







Exercise: Do example 2.17 of the text 


The idea is to set up the equations in the form : 

ax 1 + bx 2 — a 
cx 1 + dx 2 = p 

Using Cramer' s rule, the solutions are then : 
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Mechanical Rotational Systems 


• Moment of inertia: 




• Viscous friction: 



• Torsion: 






m = mn\ 





TABLE 2.5 

Torque-angular velocity, torque-ang 
rotational relationships for springs, 

-ular displacement, and impedance 
viscous dampers, and inertia 


Torque- 

Torque- 



angular 

angular 

Impedence 

Component 

velocity 

displacement 

Z M (s) = T(s)Ws) 


table 02 05 



• Model of a torsional pendulum (pendulum in clocks 
inside glass dome): 

Moment of inertia of pendulum bob denoted by J 
Friction between the bob and air by B 
Elastance of the brass suspension strip by K 



Transfer function 




Transfer Functions for Systems with Gears 




driven gear, 


Gear 2 




D 






Electro-mechanical System Transfer Functions 






Equivalence between Mechanical and 

Electrical Systems 














mass = M ► inductor 


= M hcnric 


vi scouts damper = f v »- 


spring = K 

applied force = fit) 

velocity = v(f) 




k 




voltage source = J[t) 
mesh current = v(f) 






Differential equations as mathematical models of physical 
systems: similarity between mathematical models of electrical 
circuits and models of simple mechanical systems (see model of 
an RCL circuit and model of the mass-spring-damper system). 


Ri(t ) + L 





+ B 


dx(t ) 



dt 


Kirchhofft and Newton's laws lead to mathematical models that 
describe the relationship between tie input and output of dynamic 
system. One such model is the time invariant, nth - order differenthl 
equation: 


d n c d"~ l c 

b Q i : — b ...+ Cl 

df dr 1 



+ a n c — b 

0 m 




where c(t) is the output, and r(t) is the input. 

This leads to the concept of the Transfer Function which is the corner 
stone of the control theory. Transfer functions are defined from the 
Laplace Transform. 



Summary of through- and across-variables for physical 
systems: 


System 

Variable 

Through 

Element 

Integrated 

Through 

Variable 

Variable Across 
Element 

Electrical 

Current, i 

Charge, q 

Voltage 
difference, v 21 

Mechanical 

Translational 

Force, F 

Translational 
momentum, P 

Velocity 
difference, v 21 

Mechanical 

rotational 

Torque, T 

Angular 
momentum, h 

Angular velocity 
difference, co 21 

Fluid 

Fluid volumetric 
rate of flow, Q 

Volume, V 

Pressure 
difference, P 2 , 

Thermal 

Heat flow 
rate, q 

Heat energy, 
H 

Temperature 
difference, t 21 


Integrated 
Across Variable 


Flux linkage, 

A 2 i 


Displacement 
difference, y 21 


Angular 
displacement 
difference, 0 2 


Pressure 


21 


momentum, y 21 
































I t Field current 

J constant 

1 — UUJJ — 



Armature circuit 


i 



Model of electromechanical systems. 
Model of a servomotor: 



• Model of temperature-control system: 


q e (t) = heat flow given by electric heater 
q t (t) = heat flow via liquid entering tank = 
q l (t) = heat flow into liquid 
q (t) = heat flow via liquid leaving tank = 

q (t) = heat flow (loss) via tank surface 

R e (0 + Qi (0 = Qi (0 + Qo (0 + Qs (0 



Input flow at 
temperature T f 


Liquid at 
temperature T 


Ambient temperature 
of air = T m 


Heater 


Output flow at 
temperature T 



H = specific heat of liquid 
C = thermal capacity 
R = thermal resistance 
V = volume of liquid 



+ VHT 


e 



